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ABSTRACT 

At  evening  twilight  on  4  October  1967  a  rocket  carrier  released 
three  barium  clouds  off  Wallops  Island,  Virginia.  The  vehicle  also 
carried  a  multifrequency  VHF  beacon  for  the  observation  of  radio  fre¬ 
quency  dispersive  phase  and  amplitude  variations  in  signals  propagated 
through  the  release  clouds,  and  as  a  derived  quantity,  total  electron 
content  along  the  propagation  paths.  The  first  cloud,  at  103  kilometers, 
introduced  no  measurable  dispersive  phase  variations  and  only  a  brief 
decrease  in  signal  strength.  The  other  clouds  at  187  and  226  kilometers 
produced  large  and  unusual  effects  on  both  dispersive  phase  and  signal 
strength  data. 

Analysis  of  the  dispersive  phase  data  obtained  during  rocket  descent, 
when  the  rocket  had  emerged  from  behind  the  release  clouds,  yielded  a 
normal  electron  density  profile.  A  preliminary  attempt  at  interpreting 
the  data  relevant  to  the  second  release  employed  ray  tracing  through 
spherically  symmetric,  diffusing  cloud  models  with  a  gaussian  electron 
density  distribution.  By  varying  the  parameters  defining  the  cloud,  a 
model  ’as  devised  that  closely  duplicated  both  the  dispersive  phase  and 
signal  strength  data. 
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I .  INTRODUCTION 


The  general  subject  of  and  specific  experiments  relating  to  chemical 
releases  at  very  high  altitudes  have  been  well  reviewed  in  the  literature 
(Bates,  1950;  Rosenberg  and  Golomb,  1963;  Rosenberg,  1966;  Lloyd  and 
Sheppard,  1966;  Foppl,  et  al,  1967;  Haerendel  and  Lust,  1968a;  Haerendel 
and  Lust,  1968b).  The  element,  barium,  as  a  release  constituent  has 
proven  of  particular  interest  because  it  ionizes  readily  in  sunlight  and 
the  ion  possesses  strong  resonant  scattering  lines  in  the  visible  and 
near-visible  part  of  the  spectrum.  The  formation  of  easily  observable 
barium  plasmas  at  high  altitudes  has  proven  very  useful  in  the  determina¬ 
tion  of  vector  electric  fields  present,  and  in  the  generation  of  magnetic 
field-aligned  ion  density  inhomogeneities  (tentatively  attributed  to  the 
development  of  large-scale  plasma  instabilities). 

For  experiments  in  which  a  marked  disturbance  of  the  high  altitude 

natural  environment  is  desired,  there  is  interest  in  release  particle 

8-3  3 

densities  of  perhaps  10  cm  in  a  volume  of  500  km  .  If  the  particle 
desired  is  the  barium  ion  and  an  efficiency  of  10%  is  assumed  in  convert¬ 
ing  vaporizable  barium  to  the  barium  ion  via  solar  ionization,  chemical 
releases  containing  approximately  1000  moles  (137  kg)  of  vaporizable 
barium  are  required.  A  chemical  release  of  this  magnitude  is  no  trivial 
problem.  For  diagnostic  purposes  it  becomes  important  to  know  the  release 
densities  achieved  and  the  density  distribution  in  the  release  volume. 
Appropriate  measurements  of  these  parameters  may  derive  from  observations 
of  spectral  line  (or  band)  emission  and,  for  plasmas,  from  observations 
of  radio  wave  signals  propagated  through  or  scattered  from  the  plasma 
cloud  or  from  direct  probe  measurements  in  the  plasma.  Unfortunately, 
a  barium  ion  plasma  of  the  density  and  extent  here  indicated  is  optically 
thick  for  a  substantial  period  of  time  in  the  prominent  barium  ion  tran¬ 
sitions  driven  by  solar  excitation.  For  even  the  hyperfine  emission 
arising  from  the  presence  of  barium  isotopes,  the  cloud  is  optically 
thick  at  early  times.  Thus,  radio  wave  or  direct  probe  measurements  are 
important  in  determining  the  early  time  properties  and  formation  of  plas¬ 
mas  resulting  from  chemical  releases. 
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In  late  1966,  the  Ballistic  Research  Laboratories  (BRL)  were  invited 
to  participate  in  the  chemical  release  experiment  to  be  conducted  by  the 
Air  Force  Cambridge  Research  Laboratories  (AFCRL) ,  The  experiment  en¬ 
compassed  the  launching  of  two  sounding  rockets  and  payloads,  each  con¬ 
taining  three  barium  release  canisters  for  deployment  at  different 
altitudes.  The  BRL  participation  consisted  of  instrumenting  the  nose 
section  of  each  vehicle  with  a  multi frequency  VHF  beacon  for  transmission 
downward  through  the  release  plasma,  and  for  operating  and  servicing  the 
rocket  and  ground  equipments  in  the  field.  The  experiments  were  conducted 
at  evening  twilight  on  3  and  4  October  1967,  with  vehicle  launchings 
from  Wallops  Isiand,  Virginia.  The  flights  were  given  the  code  names 
URSULA  and  VERA,  respectively.  The  chemical  releases,  six  in  number 
occurred  at  or  near  the  programmed  altitudes.  The  BRL  radiowave  experi¬ 
ment  failed  early  in  flight  on  3  October  (attributed  to  a  payload  shroud 
failure)  but  operated  successfully  during  the  flight  on  4  October.  The 
following  sections  discuss  the  nature  of  the  BRL  experiment  and  the 
results  obtained  from  the  4  October  flight  (VERA) . 

II.  THE  VHF  EXPERIMENT 

The  propagation  experiment  employed  a  three -frequency  beacon  carried 
in  the  nose  cone  of  the  chemical  release  rocket.  Significant  character¬ 
istics  of  the  transmitting  and  receiving  system  are  given  in  Table  I. 
Figures  1  and  2  show  the  rocket  payload  components  and  configuration 
similar  to  those  employed  in  the  experiment  here  described.  A  discussion 
of  very  similar  instrumentation  and  its  use  in  studies  of  the  ionosphere 
may  be  found  in  publications  of  the  BRL  (Dean,  et  al,  1968).  The  signals 
received  on  the  ground  are  electronically  processed  as  shown  in  Figure  3, 
and  the  processed  information  stored  on  paper  chart  and  magnetic  tape  for 
subsequent  analysis  and  computer  processing.  For  VERA  the  583  MHz  trans¬ 
mission  was  not  employed. 

In  operation,  the  beacon  transmits  continuously  as  the  rocket  rises. 
Upon  release  of  the  chemical,  the  transmitter  is  first  enveloped  by  the 
chemical  then  rises  above  the  release  cloud  and  transmits  through  the 
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Figure  2.  Rocketborne  Antennas  for  Four- Frequency  Beacon. 
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Figure  3.  Block  Diagram  of  Four  Frequency  Ground  Station. 
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Table  I.  Characteristics  of  Propagation  Experiment  Equipments  > 

Transmitter  I 

— * - *  ! 

Frequency:  36.44  MHz,  72„88  MHz  and  145.76  MHz  ! 

Operating  Mode:  CW  f 

Power  Output:  130  mw  at  36.44  MHz 
260  mw  at  72,88  MHz 

105  mw  at  145.76  MHz  * 

t 

Rocketbome  Antennas 

Type:  Loops  at  36.44  MHz  and  72.88  MHz 
Folded  dipole  at  145.76  MHz 

Beam  Width:  A  minimum  of  90  degrees  at  all  frequencies,  (plus 
and. .minus  45  degrees) ,  with  a  major  pattern  lobe  directed  to  the 
rear  along  the  longitudinal  axis  of  the  rocket  with  a  maximum  var¬ 
iation  of  less  than  +1.5  db  within  the  beam  width  cone. 

Polarization:  Linear  with  antenna  axis  prependicular  to  rocket  longi¬ 
tudinal  axis.  , 

l 

Gain:  (With  respect  to  a  linearly  polarized  isotropic  reference  an¬ 
tenna)  -14  db  at  36.44  MHz,  -8  db  at  72.88  MHz,  and  -4  db  at  145.76 
MHz. 

Ground  Antennas 

Tyne:  Crossed  dipoles  at  all  frequencies  at  Launch  Site;  crossed 
uipoles  at  36.44  MHz  and  72.88  MHz,  helices  at  145.76  MHz  at  Manolin 
Site. 

Polarization:  All  frequencies  had  both  right  and  left  circular  polar¬ 
ization  at  both  sites. 

Receivers 

Launch  Site:  BRL-built  combination  receivers  and  phase  locked  track¬ 
ing  loops. 

Manokin  Site:  BRL-built  receivers  and  commercial  tracking  filters. 

Resolution  (AGC) 

Launch  Site:  About  +T.5  db  at  all  frequencies. 

Manokin  Site:  About  +3.0  db  at  all  frequencies. 
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cloud  to  an  appropriately  placed  receiver  (Figure  4).  If  the  release 
cloud  is  or  becomes  ionized,  the  several  frequencies  will  experience 
differing  attenuations  and  differing  phase  path  alterations.  If  the 
radio  waves  experience  no  angular  refraction  and  the  transmissions  are 
phase  coherent  or  extremely  frequency-stable,  the  observed  phase-path 
differences  and  signal  attenuations  can  be  simply  related  to  total  elec¬ 
tron  (or  ion)  content  of  the  plasma  cloud.  In  the  more  general  case, 
aid  considering  refraction  and  finite  cloud  sizes,  the  phase-path  and 
attenuation  observations  must  be  combined  in  a  fairly  complex  computer 
program  to  furnish  both  electron  content  and  distribution,  and  then  only 
if  additional  information  on  cloud  size  and  location  is  available  to 
initiate  the  computer  routine. 


Since  a  plasma  is  dispersive  at  radio  frequencies,  the  terms  dis¬ 
persive  phase  and  dispersive  doppler  are  used  in  describing  the  real 
phase  or  frequency  difference  between  differing  frequency  signals  pro¬ 
pagated  through  such  a  medium.  Thus,  in  Figure  3,  the  quantity  DD  is  the 
Dispersive  Doppler;  it  is  a  frequency  which  is  approximately  proportional 
to  the  rate  of  change  of  integrated  electron  content  along  the  low  fre¬ 
quency  ray  path.  The  electron  content  may  vary  because  of  changes  in 
ionization  along  the  ray  path,  or  motion  of  the  rocket  thru  the  ioniza¬ 
tion,  or  a  combination  of  the  two.  If,  in  Figure  3,  the  observed  dis¬ 
persive  doppler  for  36.44  MHz  and  145.76  MHz  is  designated  $DD(fj,f2)  then, 

^Dd^I'V  =  2^1  '  (1) 


since  ^  and  f2  are  phase  coherent.  The  quantities  4^  and  <P2  are  given 


by, 


= 


2uf. 

l 


ds/ci 


(2) 


o 

where,  c^  =  propagation  velocity  of  the  ith  frequency  and  is  dependent 
upon  the  location  in  space  of  the  element  of  path  length,  ds.  In  the 
complexities  of  real  life,  the  electron  density  is  an  implicit  function 
of  the  propagation  velocity  and  the  propagation  path  depends  upon  the 
spatial  distribution  of  electron  density.  Without  some  assumptions  on 
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Figure  4.  Elevation  View  of  Rocket  Trajectory. 
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the  plasma  distribution  and/or  a  time  function  relationship  for  distri¬ 
bution  changes,  actual  plasma  electron  densities  cannot  be  derived  from 
observations  of  4^  or  <j>^.  However,  if  the  transmitter  and  receiver  are 
outside  the  plasma,  a  rather  gross  estimate  of  plasma  location  and  dis¬ 
tribution  permits  a  determination  of  total  electron  content  between 
transmitter  and  receiver  using  a  "bootstrap"  ray  trace  computer  program. 
Recognising  that  the  Appleton-Hartree  dispersion  equation  and  a  complete 
ray  trace  program  (Davies  1965,  Radio  Science,  1968)  are  required  for  a 
general  treatment,  a  simple  expression  for  c^  in  terms  of  electron  den¬ 
sity  is  possible  on  the  assumption  of  a  transmitting  frequency  large  com¬ 
pared  to  either  the  plasma  or  electron  gyromagnetic  frequencies.  Thus, 


ci  -  V  ■  c0[1- 


2 

e  N 


i-k 


me  oj 
o 


r,  80.5Ni->S  .  r,  40.3Nn-l 

ci  ■  co[1-  -pH  ■  co!l-  -pH 

l  l 


(3) 


where,  N  is  the  electron  density,  cq  is  the  vacuum  velocity  of  propaga¬ 
tion;  y  is  the  index  of  refraction  (real),  e  and  m  are  the  electron  charge 
and  mass  respectively, is  the  free  space  permittivity,  and  u>  is  the  pro¬ 
pagating  wave  angular  frequency.  Assuming,  further,  that  the  phase  paths 
at  fj  and  f2  are  identical,  substitution  of  Equations  (3)  and  (2)  into  (1) 
gives. 


,  .  80.5CK2-1)  'S 

lDB  =  ’  Kc  f. 

0  1 


Mds 


(4) 


where,  K=f2/f^.  Again  it  is  emphasized  that  in  a  real  experiment  employ¬ 
ing  two  frequencies,  the  measured  $DD  is  not  unambiguously  related  to 
Nds  as  in  Equation  (4);  however,  employment  of  additional  coherent  trans¬ 
missions  and  signal  attenuation  measurements  can  significantly  reduce  the 
ambiguity. 

The  radiowave  absorption,  per  unit  path,  is  given  by  the  expressions 
(Davies  1965,  81-82): 

Non  Deviative: 

k.  «  (e2vN)(8ir2eomc^)"1  (5) 
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Deviative: 


=  (e^vN)  [8tr^eonic£?(l-e^N/4u^meof?)Js]  *  (6) 

for  f.  much  greater  than  the  collision  (v)  and  electron  gyromagnetic  fre- 

1  _i 

quencies,  and  where  is  the  absorption  in  nepers  (8.69db)m  .  In  the 
propagation  experiments  here  discussed,  electron-neutral  and  electron-ion 
collisions  must  be  considered  in  application  of  Equations  (5)  and  (6).  At 
the  altitudes  of  interest  here,  absorption  becomes  significant  only  near 
the  plasma  frequency  and,  hence,  is  not  a  sensitive  indicator  of  electron 
density  at  other  than  the  plasma  frequency. 

As  noted  above,  interpretation  of  the  multifrequency  propagation  mea¬ 
surements  in  terms  of  electron  content/electron  density  requires  applica¬ 
tion  of  a  ray  trace  program  and  some  assumptions  on  the  location  and  elec¬ 
tron  distribution  for  the  chemical  release  plasma.  A  discussion  of  the 
ray- trace  program  is  given  in  Appendix  I;  distribution  and  location  assump¬ 
tions  are  discussed  in  later  paragraphs.  Defocussing  of  the  radio  signals 
by  the  relatively  small  release  plasma  is  significant  and  the  effect  must 
be  considered  in  interpretation  of  the  attenuation  measurements. 

III.  FLIGHT  SUMMARY 

As  indicated  earlier,  propagation  data  were  obtained  only  for  the 
second  of  the  two  rocket  probes  launched.  Trajectory  characteristics  are 
shown  in  Figures  4  and  5  and  the  release  parameters  in  Table  II  (Best, 

1968;  Noel,  1968).  Performance  of  the  vehicle,  release  locations  and 
chemical  yields  were  as  predicted.  The  propagation  experiment  functioned 
properly  with  adequate  signals  received  at  the  ground  station  located  near 
the  launcher.  Signals  at  the  remote  receiver  (Manokin  Site)  were  lost 
from  the  time  of  second  release  to  shortly  before  rocket  apogee.  The  third 
release  produced  loss  of  signal  phase  lock  for  perhaps  3  seconds  at  the 
launcher  receiver  site  and  initiated  rocket  tumble. 

IV.  DATA 

The  processed  data  from  the  propagation  experiment  are  the  dispersive 
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Figure  5.  Plan  View  of  Rocket  Trajectory. 
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TABLE  II.  Chemical  Release  Parameters,  Flight  VERA  (Oct.  4,  1967 
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phase  differences  derived  from  suitable  mixing  of  the  36.44-145.76  MHz 
signals  ($DD  of  Equation  1),  and  signal  strengths  for  the  three  propagated 
frequencies  as  measured  at  the  receiving  sites.  For  the  sake  of  brevity, 
the  beacon  frequencies  will  hereafter  be  referred  to  as  36,  73,  and  146  MHz. 
TVo  receiving  sites  were  used,  one  at  the  launch  site  and  the  other  near 
the  town  of  Manokin,  about  35  kilometers  north  and  23  kilometers  west  of 
the  launch  site.  The  dispersive  phase  (integrated  dispersive  doppler) 
measured  at  the  Manokin  site  is  shown  in  Figure  6.  At  the  time  of  re¬ 
lease  of  the  second  cloud  (133.2  seconds)  the  signal  strength  dropped  so 
that  the  tracking  filters  required  for  dispersive  doppler  measurements 
lost  lock.  The  signal  could  not  be  reacquired  until  225  seconds,  shortly 
before  peak  altitude.  Continuous  data  were  obtained  for  the  remainder  of 
the  flight. 

Dispersive  phase  from  launch  site  is  shown  in  Figure  7.  Continuous 
data  were  obtained  throughout  the  entire  flight  except  for  a  period  of 
about  three  seconds  following  the  third  releae.  The  dotted  line  is  the 
dispersive  phase  that  would  have  been  observed  had  the  rocket  been  flown 
into  an  undisturbed  ionosphere  with  the  same  profile  as  that  computed  from 
the  downleg  data  (see  below) .  The  difference  between  the  dotted  curve  and 
that  for  the  actual  data  can  therefore  be  attributed  to  the  effect  of  the 
barium  clouds.  It  is  evident  that  the  first  cloud  (71  seconds)  produced 
no  measurable  dispersive  doppler  while  the  second  and  third  events  (133.2 
and  179.0  seconds,  respectively)  produced  marked  dispersive  phase.  Effects 
due  to  the  clouds  ended  at.  about  225  seconds  and  the  dispersive  phase  due 
to  the  undisturbed  ionosphere  was  measured  for  the  remainder  of  the  flight. 
An  unusual  feature  of  the  launch  site  data  is  the  large  discontinuity  in 
dispersive  phase  which  appears  to  occur  at  the  time  of  third  release.  Such 
a  discontinuity  has  not  been  observed  prior  to  this  flight  but  it  will 
later  be  shown  that  discontinuities  can  occur  when  multipath  propagation 
is  experienced.  Further,  it  will  be  shown  that  the  short  section  of  seem¬ 
ingly  inconsistent  data  between  174  and  179  seconds  should  not  be  plotted 
as  shown.  It  should  be  replaced  by  a  discontinuity  at  174  seconds  and  be 
plotted  slightly  above  the  dotted  line  for  the  undisturbed  ionosphere. 
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Figure  6.  Manokin  Site  Dispersive  Phase  Vs  Time. 
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Figures  8  and  9  show  the  36  MHz  signal  strength  measured  using  the 
left-hand  receiver  AGC.  There  was  a  brief,  sharp  decrease  in  signal  as 
a  result  of  the  first  release  at  72  seconds,  after  which  the  signal  re¬ 
turned  to  approximately  its  undisturbed  value.  The  signal  dropout  arises 
from  a  combination  of  two  circumstances;  early  time  thermal  ionization 
in  the  thermite  burning  process  (resulting  in  absorption) ,  and  strong 
defocusing  by  the  localized  electron  cloud.  The  second  release  at  133.2 
seconds  caused  a  similar  decrease  which  was  followed  by  a  couple  of  sec¬ 
onds  of  signal  enhancement.  The  signal  enhancement  may  result  from  focus¬ 
ing  by  the  essentially  neutral  cloud  existing  immediately  after  exhaustion 
of  the  products  from  the  canisters.  The  signal  then  decreased,  with  super¬ 
imposed  oscillations,  until  174  seconds  when  it  rapidly  increased  to  an¬ 
other  period  of  signal  enhancement  until  the  third  release  at  179  seconds. 

I 

The  third  release  caused  another  brief,  sharp  decrease  in  signal  strength 
after  which  the  signal  only  partially  recovered  to  an  essentially  constant 
value,  with  superimposed  fading  for  the  remainder  of  the  flight.  Figure 
10  shows  the  signal  strength  measured  using  the  tracking  filter  AGC  for 
comparison  with  Figure  9. 

Figures  11  and  12  show  the  73  MHz  signal  strength  measured  using  the 
left-hand  receiver  AGC.  The  data  are  qualitatively  similar  to  the  36  MHz 
data  except  that  the  fading  oscillations  after  seco r'A.  release  are  more 
rapid  and  the  sudden  increase  in  signal  strength  prior  to  third  release 
occurs  earlier,  at  166.7  seconds.  The  oscillations  or  "beats"  which  are 
quite  noticeable  at  36  and  73  MHz  from  140  to  160  seconds  (from  7  to  27 
seconds  after  release)  possibly  arise  from  a  dense  core  of  ionization  or 
an  ionization  front  advancing  from  the  sunlit  side  of  the  cloud.  The 
apparent  integral  relation  of  the  beat  periods  of  the  two  frequencies 
suggests  this  interpretation. 

Figures  13  and  14  show  the  146  MHz  signal  strength  measured  using 
the  left-hand  receiver  AGC.  The  major  difference,  compared  to  the  lower 
frequencies,  occurs  between  the  second  and  third  releases.  After  the 
initial  dropout  there  is  a  gradual  decrease  and  recovery  without  the  fad¬ 
ing  that  is  characteristic  of  the  lower  frequencies.  There  is  also  no 
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sudden  increase  in  signal  strength  prior  to  third  release. 


V.  ANALYSIS 

The  computation  of  electron  density  profiles  from  dispersive  phase 
normally  assumes  either  a  horizontally  stratified,  time  invarient  iono¬ 
sphere  or  one  that  can  be  effectively  reduced  to  that  case.  For  this 
experimental  technique  the  method  for  computing  electron  density  profiles 
has  been  described  (Dean  and  Lootens,  1968).  The  effect  of  horizontal 
gradients  and  temporal  variations  has  been  discussed  (Lootens  and  Prenatt, 
1967) . 

The  VERA  flight  was  the  first  for  which  two  separated  receiving 
stations  were  employed.  As  a  result,  the  flight  provided  a  limited  test 
of  the  validity  of  the  "standard"  computational  technique.  Figure  15 
shows  the  electron  density  profiles  computed  using  the  downleg  data  from 
the  two  receiving  sites.  Also  shown  is  the  profile  computed  from  a  sweep 
frequency  ionosonde  located  near  the  launch  site  (Wright  and  Paul,  1967). 

A  qualitative  interpretation  of  a  number  of  features  of  the  disper¬ 
sive  doppler  and  signal  strength  data  obtained  during  rocket  ascent  can 
be  given  with  the  aid  of  Figure  16.  The  cloud,  with  an  index  of  refrac¬ 
tion  of  less  than  unity,  acts  as  a  diverging  lens.  Rays  passing  through 
the  cloud  are  diffused  over  a  rather  wide  area;  measured  signal  levels  at 
the  ground  would  accordingly  be  reduced.  The  prolonged  increase  in  ob¬ 
served  dispersive  phase  for  the  second  cloud  results  from  the  fact  that 
the  rays  arriving  at  the  receiver  continue  to  pass  near  the  center  of  the 
cloud  even  though  the  rocket  itself  is  emerging  from  behind  the  cloud. 
Three  regions  have  been  labeled  in  Fig.  16;  indirect,  multipath,  and 
direct.  In  the  indirect  region,  only  rays  that  have  passed  close  to  the 
center  of  the  cloud  are  received  (the  indirect  rays) .  In  the  direct  re¬ 
gion  only  rays  which  have  been  little  affected  by  the  cloud  are  received 
(the  direct  rays).  In  the  multipath  region  both  direct  and  indirect  rays 
reach  the  receiver.  As  the  rocket  moves  upward  and  to  the  right,  the  in¬ 
side  edge  of  the  multipath  region  moves  to  the  left  so  that,  if  a  multi- 
path  region  exists  at  all,  there  occurs  a  time  for  a  given  location  when 

34 


% 

■J 


INDIRECT" 

MULTIPATH 

"DIRECT 

RAYS 

RAYS 

ONLY 

ONLY 

Figure  16.  Sketch  Illustrating  The  Effect  of  an  Ionized. 
Cloud  on  Propagation  Paths. 


the  direct  signal  is  first  received.  At  that  time  for  a  receiver  posi¬ 
tion  corresponding  to  "Location  A"  in  Figure  16,  some  focusing  of  the 
rays  occurs  at  the  inner  edge  of  the  multipath  region  leading  to  signal 
enhancement.  The  "earliest  multipath"  effect  can  explain  both  the  sudden 
increases  in  signal  strength  at  36  and  73  MHz  and  the  sudden  change  in 
slope  of  the  dispersive  phase  at  174  seconds.  At  that  itme,  the  36  MHz 
phase-locked  tracking  filters  lost  the  weak  indirect  signal  and  locked 
onto  the  much  stronger  direct  signal  (labeled  "Ray  Path  A"  in  Figure  16). 
From  Figure  16  it  is  evident  that  between  174  and  179  seconds  the  mea¬ 
sured  dispersive  phase  is  along  the  direct  ray  path  which  is  relatively 
little  affected  by  the  cloud.  As  a  result,  the  dispersive  phase  during 
that  period  should  not  be  plotted  as  shown  in  Figure  7,  but  should  appear 
slightly  above,  and  approaching,  the  dotted  line  corresponding  to  the 
undisturbed  ionosphere. 

Following  the  third  release,  the  absence  of  effects  characteristic 
of  "earliest  multipath"  indicates  that  multipath  was  not  observed  for 
that  cloud. 

A  preliminary  quantitative  interpretation  of  the  dispersive  doppler 
and  signal  strength  data  has  been  attempted  while  the  signal  was  being 
affected  by  the  second  barium  cloud.  A  spherically  sysmmetric  cloud 
model,  whose  defining  parameters  are  shown  in  Figure  17,  was  postulated. 

The  electron  density  wa s  assumed  to  be  gaussian  with  a  total  electron 
inventor)''  beginning  at  zero  and  asymptotically  approaching  some  limiting 
value.  The  square  of  the  gaussian  scale  radius  was  assumed  to  increase 
linearly  from  some  initial  value  as  the  cloud  grew  by  diffusion.  Ion  cloud 
drift  was  determined  from  optical  data  (Noel,  1968).  A  ray  tracing  pro¬ 
gram  (see  Appendix  I)  was  employed  to  compute  the  dispersive  phase  cor¬ 
responding  to  the  cloud  models.  The  parameters  defining  the  cloud  were 
adjusted  until  agreement  was  achieved  between  computed  dispersive  phase 
and  that  actually  measured. 

Two  approaches  were  used  adjusting  cloud  parameters .  One  involved 
adjusting  the  ionization  time  constant,  x,  the  total  inventory  limit 
TIji  ,  the  initial  scale  radius,  h  ,  and  the  cloud  growth  rate,  ^^(=4D) , 
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to  bring  the  computed  and  measured  dispersive  phase  into  agreement.  Using 
this  approach,  three  clouds  were  defined,  each  of  which  reproduced  the  mea¬ 
sured  dispersive  phase  to  within  1-2  percent  between  140  and  174  seconds. 
The  several  cloud  models,  in  terms  of  scale  radius  and  total  ion  inventory, 
are  shown  in  Figure  18,  where  the  crosses  connected  by  dashed  lines  indi¬ 
cate  cloud  parameters  at  several  times  following  release  (release  occurred 
at  133.2  seconds).  The  defining  cloud  parameters  are  listed  in  Table  III. 

The  other  approach  in  determining  cloud  parameters  considered  only 
scale  radius  and  total  inventory.  A  value  of  scale  radius  was  selected 
and  the  values  of  total  inventory  that  would  reproduce  the  data  were  com¬ 
puted  for  several  times.  This  approach  resulted,  for  cloud  II,  in  the 
family  of  nearly  straight  lines,  each  for  a  different  time,  shown  in  Fig¬ 
ure  18. 

Based  on  dispersive  phase  data  alone,  it  was  not  possible  to  decide 
among  the  various  cloud  models,  each  of  which  reproduced  the  measured 
data.  To  resolve  the  ambiguity  it  was  necessary  to  include  analysis  of 
the  signal  strength  data.  As  discussed  earlier,  the  large  increases  in 
the  strength  of  the  36  and  73  MHz  signals  prior  to  third  release  were 
assumed  to  correspond  to  first  arrival  of  the  direct  radio  wave  signal. 

Each  of  the  clouds,  defined  in  Table  III,  was  considered  at  174  seconds 
and  families  of  36  MHz  rays  were  traced  through  each  cloud.  Plots  of  these 
ray  traces  are  shown  in  Figure  19.  The  family  of  rays  for  the  large  cloud 
shows  that,  at  174  seconds,  the  inside  edge  of  the  multipath  region  had 
not  yet  reached  the  observer  so  the  increase  in  signal  strength  would  not 
have  occurred  as  early  as  174  seconds.  TV'  family  of  rays  for  the  small 
cloud  shows  that  multipath  was  already  being  experienced  prior  to  174 
seconds.  The  family  of  rays  for  the  medium  cloud  shows  that  the  inside 
edge  of  the  multipath  is  almost  coincident  with  the  observer. 

The  same  exercise  was  carried  out  for  the  73  MHz  sig  il  at  166.7 
seconds.  The  results,  shown  in  Figure  20,  are  essentially  the  same  as 
for  36  MHz  except  that  for  the  large  cloud  multipath  is  not  even  possible. 
As  a  result,  it  appears  that  the  medium  cloud  exhibits  the  best  agree¬ 
ment  with  the  measured  data. 
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Figure  19  Raysets  at  36.44  MHz  for  Three  Model  Clouds  at  174.0  Seconds 
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TABLE  III.  Parameters  Defining  Cloud  II  Models  That  Reproduced  Measured  Data 
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By  a  technique  briefly  outlined  in  Appendix  I,  reductions  of  signal 
levels  by  the  defocusing  mechanism  were  computed  at  selected  times  for 
the  several  cloud  models.  Results  of  the  computation  are  given  in  Table 
IV  and,  for  the  36  MHz  ordinary  mode  signals,  are  plotted  in  Figure  10. 
Comparison  of  the  computed  and  measured  signal  strengths  at  36  MHz  would 
seem  to  indicate  a  cloud  somewhat  smaller  than  cloud  model  B  and  depart¬ 
ing  from  spherical  symmetry  after  160  seconds. 

An  attempt  was  made  to  perform  an  analysis  of  the  Cloud  III  data 
similar  to  that  done  for  Cloud  II.  However,  at  the  time  of  this  writing, 
no  satisfactory  combination  of  parameters  has  been  found  for  a  sperically 
symmetric  cloud  that  will  reproduce  the  cloud.  It  appears  that  analysis 
of  the  Cloud  III  data  will  require  use  of  a  more  elaborate  model,  for 
example  a  field-aligned,  prolate  spheriod.  Current  plans  call  for  analy¬ 
sis  of  the  data  from  both  clouds  using  such  a  model. 

VI.  CONCLUSIONS 

The  preceeding  sections  describe  a  specific  barium  release  ion  cloud 
in  terms  of  peak  electron  density,  total  ion  inventory,  ion  diffusion  rate 
and  ion  cloud  radius  for  the  interval  i0<t<40  seconds  following  release. 
The  model  for  the  cloud  assumes  a  gaussian  distribution  for  barium  ions 
and  diffusion  radially  symmetric  in  the  time  interval  examined.  The  model 
parameters  are  forced  to  fit  dispersive  phase  and  signal  strength  data 
derived  from  a  radio  wave  technique  employing  a  multi frequency  phase-co¬ 
herent  beacon  transmitting  through  the  ion  cloud  as  it  forms  and  grows. 

For  the  release  mixture  given  in  Table  I  and  deployed  at  187  km.,  the 
release  cloud  characteristics  may  be  summarized  as  follows: 

Time  Constant  for  Ba  II  formation  —  15  sec 

Diffusion  Coefficient  for  Ba  II  (approx.  perpendicular  to  B)  -- 

7.SxlO~3km23ec”1 

6  -3 

Peak  Electron  Density  (at  release+38  sec)  ----9x10  el  cm 

23 

Total  Ion  Inventory  (release+^J  -  7.0x10 

Scale  Radius  Perpendicular  to  $  (release+40  sec)  -  2.3S  km 
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TABLE  IV.  Computed  Signal  Loss  by  the  Defocusing  Mechanism  (36  MHz,  Ordinary  Mode) 
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The  derived  parameters  of  the  barium  cloud  appear  quite  reasonable,  al¬ 
though  comparative  optical  data  arc  very  sparse  at  the  early  times  at 
which  the  radio  frequency  measurements  are  made.  The  total  ion  inventory 
derived  represents  an  ion  yield  of  2.6%  of  the  total  chemical  payload 
mass.  The  diffusion  coefficient  appears  low  (approximately  one-fifth  the 
neutral  barium  atom  diffusion  coefficient  at  the  release  altitude)  but 
is  not  unreasonable,  since  the  coefficient  measured  represents  diffusion 
normal  to  the  magnetic  field  because  of  the  geometry  of  release  and  the 
radio  wave  ray  paths  involved. 

The  technique  of  utilizing  radio  wave  transmissions  to  describe 
barium  cloud  (or  any  plasma)  properties  has  several  domains  of  applica¬ 
tion.  If  ionization  time  constant  and  time  for  neutral  cloud  stabiliza¬ 
tion  are  short  compared  to  the  time  the  receiver  site  is  eclipsed  by  the 
cloud,  the  technique  is  useful  in  describing  initial  cloud  development, 
ion  cloud  diffusion,  ion  inventory,  and  peak  cloud  electron  densities. 

If  the  transmissions  through  the  cloud  take  place  at  late  time  (i.e., 
after  the  plasma  cloud  striates),  the  technique  furnishes  information  on 
cloud  structure,  signal  defocusing  and  on  electron  densities  if  a  reason¬ 
able  cloud  model  can  be  postulated. 

Perhaps  the  most  interesting  result  of  the  specific  experiment  dis¬ 
cussed  here  is  the  strong  defocusing  of  the  radio  signals  as  evidenced 
in  Figures  10,  12,  and  14.  The  signal  attenuation  is  entirely  the  re¬ 
sult  of  the  defocusing  mechanism;  absorption  plays  no  part  at  the  altitude 
of  release.  Under  circumstances  for  which  the  transmitter  is  far  removed 
from  the  plasma  cloud,  the  affected  area  on  the  ground  would  be  restricted 
to  the  size  of  the  eclipsing  plasma.  For  transmissions  from  the  ground, 
defocusing  effects  would  be  widespread  for  transmitter-receiver  distances 
large  compared  to  the  transmitter-cloud  distance. 

Considerably  more  analysis  is  required  to  understand  the  details 
(periodic  fading)  in  the  signal  strength  data. 
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APPENDIX  I:  RESUME  OF  THE  RAY  TRACE  PROGRAM 


The  ray  tracing  program  used  here  is  a  combination  of  two  earlier 
programs.  The  actual  stepping  of  the  ray  is  accomplished  by  numerical 
integration  of  a  set  of  differential  equations  similar  to  Hazelgrove's 
equations  (Hazelgrove,  1954 ;  Jones,  1966).  The  equations  are: 
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where  r,  0,  and  <p  are  the  spherical  coordinates  of  a  point  on  the  ray 
path,  Vr,  Vq,  and  are  the  components  of  the  wave  normal  direction, 
and  p  and  p'  are  the  phase  and  group  refractive  indices.  Integration  of 
the  first  three  equations  determines  the  location  of  the  ray  path  in  space 
while  integration  of  the  last  three  equations  determines  the  components 
of  the  wave  normal  direction.  Additional  equations  enable  calculation  of 
phase  path,  absorption,  and  other  parameters  which  may  be  integrated 
along  the  ray  path.  The  independent  variable,  t,  is  group  path  rather 
than  phase  path  as  used  by  Hazelgrove. 


The  remainder  of  the  ray  trace  program,  i.e.,  the  initial  aiming 
of  the  ray,  the  homing  of  the  ray  on  the  observer,  the  computation  of 
the  dispersive  phase,  etc.,  employs  some  of  the  techniques  from  earlier, 
Snell's  Law  ray  tracing  programs  (R.S.  Lawrence  and  D.  J.  Posakony,  1961; 
G.  A.  Dulk  and  W.  A.  Dean,  1962).  When  tracing  the  so-called  "indirect 
rays",  the  initial  aiming  direction  is  toward  a  point  slightly  below  the 
cloud  center.  For  "direct  rays"  the  initial  direction  is  toward  the  ob¬ 
server.  In  each  case,  the  ray  is  traced  until  it  reaches  the  "impact 


49 


plane",  a  plane  through  the  observer  and  perpendicular  to  the  transmitter- 
observer  line.  The  distance  by  which  the  impact  point  misses  the  observer 
(the  miss  distance)  is  used  to  compute  a  new  aiming  direction  and  another 
ray  is  traced.  This  process,  called  homing,  is  repeated  until  the  miss 
distance  is  less  than  a  predetermined  amount. 

Dispersive  phase  is  computed  by  averaging  the  ordinary  and  extra¬ 
ordinary  phase  paths  at  the  36  and  the  146  MHz  frequencies.  The  phase 
path  average  at  36  MHz  is  multiplied  by  four  and  subtracted  from  the 
146  MHz  phase  path  average.  The  result  is  multiplied  by  two  and  this 
quantity  is  termed  dispersive  phase. 

The  changes  in  signal  strength  as  a  result  of  focusing  and  de- 
focusing  are  measured  in  the  following  way.  A  ray  at  the  frequency 
under  consideration  is  traced  through  the  unperturbed  ionosphere  until 
it  homes  on  the  observer.  The  starting  direction  of  the  ray  is  then 
deviated  by  1/10  milliradian  in  each  of  three  mutually  perpendicular 
planes.  The  three  rays  so  defined  are  traced  to  the  impact  plane  and 
the  area  of  the  triangle  whose  vertices  are  the  impact  points  is  com¬ 
puted.  The  entire  process  is  then  repeated  with  the  electron  density 
perturbation  due  to  the  barium  cloud  included.  The  ratio  of  the  triangle 
areas  with  and  without  the  cloud  is  computed  and  expressed  in  decibels 
as  a  measure  of  signal  reduction  or  enhancement  due  to  the  cloud. 
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